Vertex-centred discretization of multiphase compositional Darcy flows on general meshes

This paper concerns the discretization on general 3D meshes of multiphase compositional Darcy flows in heterogeneous anisotropic porous media. Extending Coats’ formulation [15] to an arbitrary number of phases, the model accounts for the coupling of the mass balance of each component with the pore volume conservation and the thermodynamical equilibrium, and dynamically manages phase appearance and disappearance. The spatial discretization of the multiphase compositional Darcy flows is based on a generalization of the Vertex Approximate Gradient scheme (VAG), already introduced for single phase diffusive problems in [24]. It leads to an unconditionally coercive scheme for arbitrary meshes and permeability tensors. The stencil of this vertex-centred scheme typically comprises 27 points on topologically Cartesian meshes, and the number of unknowns on tetrahedral meshes is considerably reduced, compared with usual cell-centred approaches. The efficiency of our approach is exhibited on the nearwell injection of miscible CO2 in a saline aquifer taking into account the vaporization of H2O in the gas phase as well as the precipitation of salt. R. Eymard University of Paris-Est, France, E-mail: [email protected] C. Guichard University of Nice Sophia Antipolis, France, E-mail: [email protected] R. Herbin University of Aix-Marseille, France, E-mail: [email protected] R. Masson University of Nice Sophia Antipolis, France, E-mail: [email protected]